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than the ones obtained from integration-by-parts. 
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1. Introduction 

The mathematical structures behind Feynman integrals are fascinating and by far not yet com- 
pletely understood. Feynman integrals evaluate to transcendental functions and a straightforward 
question is which classes of transcendental functions occur in the computation of Feynman inte- 
grals. If we restrict ourselves to one-loop integrals, the answer is known: We encounter only two 
transcendental functions, the logarithm and the dilogarithm which are defined by 



U\(x) = -ln(l -x) = £ — , 

n— 1 ^ 

OO ft 

Li 2 (*) = £^, (1.1) 

Jm^ 1-1- 



and whose arguments are algebraic functions of the momenta and the masses. The logarithm and 
the dilogarithm have generalisations, the most obvious one is given by the polylogarithms defined 
by 



n=\ 



At the next stage of generalisation one encounters multiple polylogarithms defined by [jj], ^] 



l \ ^ — ' 1 2 k o 

Li»ii,m2,...,mt(-''l>-''2) ■■•i^k) = 7, mi tm ' mt- (1-3) 



Bl>«2>...>/Jjfc>0 n i 



The class of multiple polylogarithms defined by eq. ( |1.3| ) plays an important role in the calculation 
of Feynman integrals beyond one-loop. Indeed, many of the known two-loop amplitudes with 
massless particles can be expressed in terms of these functions. The multiple polylogarithms have 
a rich algebraic structure: There are two Hopf algebras, one with a shuffle multiplication induced 
from the integral representation and the other one with a quasi-shuffle multiplication induced from 
the sum representation. In addition there are convolution and conjugation operations [|[]. 

However, it is also known that there are integrals which cannot be expressed in terms of mul- 
tiple polylogarithms. The simplest example is given by the two-loop sunset integral with three 
internal non-zero masses. The corresponding Feynman diagram is shown in fig. ([I]). It is therefore 
worth to study this integral in order to learn more about the functions beyond multiple polylog- 
arithms associated to Feynman integrals. The two-loop sunrise integral has received in the past 



significant attention in the literature [|4|-|ljj. Despite this effort, an analytical answer in the general 
case of unequal masses is not yet known. In the special case where all three internal masses are 
equal, a second-order differential equation in the external momentum squared and its analytical 
solution are known [jl^. The analytical solution for the equal mass case involves elliptic functions. 
In the general case of unequal masses integration-by-parts identities [ |l7[ |T8| ] can be used to relate 
integrals with different powers of the propagators. In the case of the sunrise topology with unequal 
masses all integrals can be expressed in terms of four master integrals plus simpler integrals. This 
results in a coupled system of four first-order differential equations for the four master integrals 
[^]. This is however not yet the simplest form for the differential equations governing the two-loop 
sunrise integral. In a recent publication [ p^ ] we showed - using methods of algebraic geometry - 
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Figure 1: The two-loop sunrise graph. 

that also in the unequal mass case there is a single second-order differential equation. In this talk 
we review the derivation of the second-order differential equation. 

Algebraic geometry studies the zero sets of polynomials. A simple example is given by the 
equation 

X\X2 +*2*3 + *3*1 =0. (1.4) 

The zero set of this equation defines an algebraic variety. In this example it is easily observed 
that given any solution (x[' ,j4 ) also the point (Xxf\ Xx^\ Xx^ ) is a solution. Therefore 
eq. (L4) defines an algebraic variety in the projective space P 2 . We may then study integrals of the 
form 

/ A5(l-t^3 — ^— , (1.5) 

Xi>0 v '~ l > 

where the polynomial defining the algebraic variety appears in the denominator. This leads imme- 
diately to the question what happens in the points (1,0,0), (0,1,0) or (0,0,1), which lie on the 
boundary of the integration region and where the polynomial in the denominator vanishes. The in- 
terplay between the algebraic variety defined by the zero set of the denominator with the integration 
boundary will be the Leitmotiv in our treatment of the two-loop sunrise integral. 

2. The two-loop sunrise integral 

The two-loop sunrise integral is given in D-dimensional Minkowski space by 

2\ / 2\3-o fd D kid D k 2 1 ,„ 1S 
S(D,p 2 ) = (jU 2 ) / — ^ — f- -, r-. (2.1) 



Here we suppressed on the l.h.s. the dependence on the internal masses m\,mi and and on the 
arbitrary scale \l. It is convenient to denote the momentum squared by t = p 2 . We can trade the 
integration over the loop momenta for an integration over Feynman parameters. The integral over 
the Feynman parameters depends then on two graph polynomials [ |20| ] and reads 



S(D,t) =r(3-D)( A t 2 ) 3 - D j -^jyCO, (2.2) 
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where the two Feynman graph polynomials are given by 

& = — x\X2XT,t + (xim 2 +X2W2 + X3ra 2 ) % , 9/ = x\X2 +X2X3 +X3X1. (2.3) 
The differential two-form ft) is given by 

ft) = x\dx2 Adxi, + A dx\ +x$dxi Adx2- (2.4) 

The integration is over 

a = {[xi :x 2 :jc 3 ] € P 2 |x,- > 0,/ = 1,2,3} . (2.5) 

It is simpler to consider this integral first in D = 2 dimensions and to obtain the result in D = 4 — 2e 
dimensions with the help of dimensional recurrence relations [^l], 22]. In two dimensions this 
integral is finite and given by 

2 f CO 



5(2,0 = tfj (2.6) 

a 

In two dimensions the integral depends only on the second Symanzik polynomial & . Note that 



eq. ( |2.6| ) is similar to eq. ( |1.5| ). 

Now let us turn to Hodge structures. Hodge structures have their origin in the study of compact 
Kahler manifolds. There one has the following decomposition of the cohomology groups 

H k (X,C)= H p - c '(X), HP'9(X)=H q ' p (X). (2.7) 

p+q=k 

For a fixed k this provides an example of a pure Hodge structure of weight k. Algebraic varieties 



have a generalisation of this structure, which are called mixed Hodge structures [|23J, [24J]. In addi- 
tion, we can consider a family of Hodge structures, parametrised by a manifold [^, pq ]. This is 
called a variation of a Hodge structure. Suppose that the cohomology groups are finite dimensional. 
It follows that if a Hodge structure varies smoothly with some parameters, then 

dim# M (2.8) 

remains constant. In the following we will relate the order of the differential equation for the two- 
loop sunrise integral to the dimension of a cohomology group. If the variation with the internal 
masses is smooth and if the integral has a second-order differential equation in the equal mass case, 
it follows that there must be also a second-order differential equation in the unequal mass case. 



Now let us come back to the integral in eq. (|2.6|). From the point of view of algebraic geometry 



there are two objects of interest: On the one hand the domain of integration o and on the other hand 
the algebraic variety X defined by the zero set of & = 0. The two objects X and a intersect at the 
three points [1 : : 0], [0 : 1 : 0] and [0:0: 1]. This is shown in fig. (§). We blow-up P 2 in these 
three points and we denote the blow-up by P. We further denote the strict transform of X by Y and 
the total transform of the set {X1X2X3 = 0} by B. With these notations we can now consider the 



mixed Hodge structure (or the motive) given by the relative cohomology group [ 27 ] 



H 2 (P\Y,B\BnY). (2.9) 
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Figure 2: The intersection of the domain of in- 
tegration a with the zero set X of the second 
Symanzik polynomial. 




Figure 3: The elliptic curve y 2 = x 3 — x + 1 . 



In the case of the two-loop sunrise integral considered here essential information on H 2 (P\Y,B\BPi 
Y) is already given by H l (X). We recall that the algebraic variety X is defined by the second 
Symanzik polynomial: 



-Xix 2 x 3 t + [x\m\ + x 2 m\ + X3WI3) (x\x 2 + x 2 x 3 + X3X1) = 0. 



(2.10) 



This defines for generic values of the parameters t, m\, m 2 and 013 an elliptic curve. The elliptic 
curve varies smoothly with the parameters t, m\, m 2 and 7713. By a birational change of coordinates 
this equation can brought into the Weierstrass normal form 



y 2 z-x 3 -a 2 {t)xz 2 -a 3 {t)z 3 = 0. 



(2.11) 



The dependence of a 2 and 03 on the masses is not written explicitly. In the chart z = 1 this reduces 
to 



y 2 -x 3 -a 2 (t)x-a 3 (t) 



0. 



(2.12) 



The curve varies with the parameter t. An example of an elliptic curve is shown in fig. (|3|). It is 
well-known that in the coordinates of eq. ( 2.12 ) the cohomology group H 1 (X) is generated by 



dx , . d 
77 = — and 77 = — 77. 

y 



dt 



(2.13) 



Since H (X) is two-dimensional it follows that f) = ^r 7 ? must be a linear combination of 77 and 77. 
In other words we must have a relation of the form 

Po(t)i)+Pi(t)fl+P2(t)ri = 0. 
The coefficients po(t), p\{t) and p 2 (t) define the Picard-Fuchs operator 



(2.14) 



l(2) = P0 ^h +Pl ^di +P2 ^ 



(2.15) 
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Applying the Picard-Fuchs operator to our integrand gives an exact form: 

(0 \ 



dp. 



The integration over a yields 



L®S(2,t) = At 2 J dp = \i 2 J p 
<y da 



(2.16) 



(2.17) 



The integration of /3 over da is elementary and we arrive at 

5(2,0 



M MO- 



(2.18) 



This is the sought-after second-order differential equation. The polynomials pj(t) are given by 



Po(t) 
MO 



(mi + m 2 + WJ3J 



t — (m\ — ni2 + m^y 



MO 

MO 



t — (—mi + m2 + mj) 
t-(mi+m 2 -m 3 ) 2 \ [3t 2 - 2M {m t - M 20 q + 2M y 10 ] , 
9t 6 - 32M 100 ? 5 + (37M 200 + 70Mu ) t 4 - (8M 300 + 56M 2 i + 144M m ) t 3 

- (13M 400 - 36M 3 io + 46M 220 - 124M 2 n) ? 2 

- (-8M500 + 24M 4 io - I6M320 - 96M 31 1 + 144M 22 i ) f 

- (M 600 - 6M 510 + 15M 420 - 2OM330 + I8M411 - 12M 321 - 6M 222 ) , 
3t 5 - lM m t A + (2M 200 + 16M 110 ) t 3 + (6M300 - 14M 210 ) t 2 

- (5M 400 - 8M310 + 6M220 - 8M211 ) t + (M500 - 3M 410 + 2M 320 + 8M31 1 - IOM221 ) , 

- 18? 4 + 24M im t 3 + (4M 200 - 40Mi 10 ) t 2 + (-8M300 + 8M 2 i + 48M! n )t 

2 

IT 

+ (-2M 4 oo + 8M 3 io- 



11% 

12M 22 o - 8M211 ) + 2c (t,mi , m 2 ,m 3 )ln -± 



2 2 

m 2 m 3 

+2c(f,m2,m3,mi)ln— j + 2c (t,m 3 ,mi,m2)ln- ~ 



(2.19) 



with 



c(t,m l ,m 2 ,m 3 ) = 

(-2m 2 + m\ + m 2 ) t 3 + (6m\ - 3m 2 - - 7m 2 m 2 - lm\m\ + 14m 2 m 2 ) t 2 

/ 6 6 6 4 2 4 2 24 24 42 2 4\ 

+ (—6m! + 3m 2 + 3m 3 + \ \mim 2 + llm 1 m 3 — 8m]m 2 — SmjT^ — 3m 2 m 3 — 3m 2 m 3 J ? 

+ (2mi - m 2 - - Smfm 2 - 5mf m 2 + m 2 m 2 + m 2 m3 + 4m 2 m 2 + 4m 2 m3 

+3m jm 2 + 3m\m\ — 6m 2 m 3 + 2m\m\m\ — m\m\m\ — m 2 m 2 m 3 ) . (2.20) 

In order to present the result in a compact form we have introduced the monomial symmetric 
polynomials in the variables m 2 , m^ and m 2 . These are denned by 



M 



A] A2 A3 



m, 



1m 



(2.21) 



where the sum is over all distinct permutations of (^1,^2,^3). In the equal mass case eq. (2.18) 
reduces to the well-known result of [13]. The differential equation in eq. ( 2.18j ) has been confirmed 
with numerical methods in [|28t] . 
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